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Abstract 

A numerical renormalization group technique recently developed by one of 
us is used to analyse the Coulomb pseudopotential (/i*) in Cgo for a variety 



^ ■ of bare potentials. We find a large reduction in fi* due to intraball screening 

^r ■ alone, leading to an interesting non-monotonic dependence of fi* on the bare 

interaction strength. We find that /x* is positive for physically reasonable 
0^ . bare parameters, but small enough to make the electron-phonon coupling a 



viable mechanism for superconductivity in alkali-doped fullerides. We end 



i-rt ■ with some open problems. 
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The discovery of the icosahedrally symmetric molecule Ceo in 1985 ^, and the subsequent 
discovery of the many interesting electronic properties of doped fuUerene lattices has led to 
an explosion of research work on these compounds. In this Letter we will be concerned with 
the effects of Coulomb interactions on fullerene superconductivity. 

Cqo forms an FCC lattice which can be doped with up to 6 alkali metal atoms per 
Ceo- The triply-doped compounds A^Cqq (where A = K, Rb, Cs) are bad metals but 
superconduct at relatively high temperatures (~ 20K for i^aCeo) HFH- They also have a very 
short coherence length of ^ ~ 30A, which is only about 2 lattice spacings of the FCC lattice. 
Experimental evidence ||5|-|10| , supported by numerous phonon calculations ||TT|-[18| , indicates 
that some of the intraball phonon modes are very high in energy and sufficiently strongly 
coupled to the electrons to account for these high T^ s. Most of the phonon calculations 
assume that the Coulomb repulsion between the electrons comprising a pair are small, as 
they are in a usual metallic superconductor. (Effective Coulomb interactions are typically 
parameterized by a dimensionless number /i* = N{Ep)U* |[T^, where N{Ep) is the density 
of states at the Fermi Surface (FS) and U* is an effective interaction matrix element between 
electrons at the FS. For a metal /i* ~ 0.2). However, there are two significant differences in 
the case of superconducting fuUerides which make this correspondence uncertain. Firstly, 
the coherence length, which is a measure of how close the electrons in a pair are, is an order 
of magnitude shorter in AsCeo than in metallic superconductors. Furthermore, in metallic 
superconductors the effective Coulomb interaction between the electrons comprising a pair 
is reduced by the highly retarded nature of the phonon-mediated attraction, as shown by 



Anderson and Morel in 1962 19 



where ujd is the Debye frequency, W is the electronic bandwidth, and /i measures the instan- 
taneous Coulomb repulsion. In metals Ud is typically two orders of magnitude lower than 
W , which accounts for the small /i*. However, the electronic bandwith of the conduction 
level of A3C60 is about W ~ 0.5el^, while the phonon modes go up to an energy of Q.2eV 



0], which makes retardation highly ineffective in reducing /x*. 

However, it has been proposed by Baskaran and Tosatti |^, and simultaneously by 
Chakravarty and Kivelson [^, that a repulsive bare interaction between the electrons on 
a single fullerene molecule can result in an effective attraction between electrons in the 
conduction level. Second-order perturbation theory (PT) in U /t in the tight-binding model 
on the truncated icosahedral (TI) lattice with the on-site Hubbard interaction U shows 



that jj* becomes negative beyond U/t ~ 3 |^T|. Exact diagonalizations on small clusters 



2^ , |23[| also demonstrate the possibility of /x* < for intermediate U/t. However, increasing 



the range of the interaction increases /x* in second-order PT [P4| , P5| , and the validity of 
second-order PT is questionable in the region where /i* < p5 |. 



A numerical RG method has been developed recently by one of us to address these issues 
and get a reliable estimate of /i* p6|. This is based on the extention of Wilson's momentum 



shell integration |^ to the case of condensed fermion systems as developed by Shankar 
P8|| . The crucial distinction between the RG as applied to critical phenomena in statistical 
mechanics, and the RG as applied to condensed fermion systems is that in the latter case 
there are an infinite number of flowing coupling constants which are relevant to the low- 
energy behavior of the system (in infinite volume) [^ . To illustrate the method for finite 
sizes, let us consider the case of interacting electrons on Ceo- We start with the Hamiltonian 

H = -tY. (49s + h.c.) + Y. V{i,j)n,nj (2) 

{i3,s) hi 

where the first term is just the tight-binding approximation and the second term represents 
a generic density-density interaction between the electrons. (We have also performed calcu- 
lations for two different hopping matrix elements for the purely pentagonal versus the other 
bonds [Q, but since the results are almost identical to the ones for the above hamiltonian, 
we report only the latter here). We solve the tight-binding problem and get 60 eigenstates, 
each belonging to a representation of the icosahedral group Ih [^ • We can then go to the 
zero-temperature fermionic path integral with the action 



/oo 
~ Zl / 77r-T^'^(^l+^2-^3-tU4)V'({fcj})^fcis%2s'%33"7fc4S (3) 

where 77, r/ are Grassmann variables representing the fermions, and Ks'(^i5 ^2, ^3, ^4) is the 
matrix element of the bare interaction, and Z is the quasiparticle residue, or the wave 
function renormalization. For the bare theory Z = 1 and the energies are the tight-binding 
energies. This method produces results in excellent agreement with exact diagonalizations 



for small clusters 26 



We now integrate out the variables corresponding to the high energy levels (levels far 
from the conduction level) step by step. We go to two-loop order in the interaction matrix 
elements to compute effective values for the energy, the wave function renormalization, and 



the interaction. Each integration generates effective interactions for the remaining levels [26 



We continue until only the conduction level remains, and obtain an effective hamiltonian 



for the conduction level. In accordance with the formalism for infinite fermion systems |^ 
we allow all the independent couplings to flow separately. In order to reduce computation 
time we use the full symmetries of Ih to keep only symmetry-reduced matrix elements, in 
analogy with the Wigner-Eckhart theorem. This still leads to ~ 10^ flowing couplings. We 
also ignore the frequency dependence of the effective interactions (they are irrelevant in the 
RG sense for properties near the fermi surface |]28|). Apart from keeping track of Z, going 
to two-loop order enables us to generate an internal validity criterion. As far as the effective 
interactions in the conduction level are concerned, the smallest energy denominators come 
from the excitations to levels close to the conduction level, which will be integrated out 
during the final steps of the RG. We compare the two- loop to the one- loop contribution to 
the effective interactions at the final RG step. When their ratio becomes ^ 1 we know that 
all loops will be important, and that truncating to two-loop is invalid. As will be evident 
from our numerical results, this method is reliable up to intermediate values of U /t ~ 5, 
which fortunately includes the physical regime for Cgo §]• 

We have considered two classes of models for the bare potential; the extended Hubbard 



model with on-site and nearest neighbor interactions, a screened Ohno potentiaL We also 
considered the problem of only two additional electrons in the Ti^ conduction level (as long as 
the three- and higher-body interactions are negligible compared to the two-body interaction, 
these results should hold for arbitrary numbers of electrons). The two additional electrons 
can be in one of three channels Ag, Tig, Hg, which roughly correspond to total angular 
momentum L = 0, 1, 2 respectively |^. Let us now proceed to the results. 



In Fig. m we show the effective interaction U* /t for the on-site Hubbard model as a 
function of the bare interaction U/t. From the non-monotonicity it is very evident that no 
naive correspondence can be made with the Anderson- Morel formula, eq(l). For U/t < 5 
the triplet Tig pair state has the lowest energy. When considered in the context of a lattice 
problem, this will lead to itinerant ferromagnetism. There is also a level crossing between 
the Ag and the Hg pair states at U/t ~ 3.5, and for U/t > 5 the Ag state is the ground 
state. Fig. ^ shows how far one can trust these results by plotting the ratio of the third 
order to the second order at the final step of RG as a function of U/t. We see that these 
results are reliable up to U/t ^ 5 — Q. For U/t > 6 the Ag channel shows negative U*, 
which was precisely the original counterintuitive claim that inspired this work. However, it 
unfortunately happens in a region where we cannot trust the result. Fig. ^ compares various 
different approximations for U*/t for the Ag channel, from which it is clear that perturbation 
theory fails for physically interesting U/t ^3 — 5, and even one- loop RG differs significantly 
from the two-loop result. Finally, the ladder resummation |3^, which corresponds in our 
language to a one-loop RG with only the Cooper channel, is also significantly different from 
the two- loop result. Qualitatively similar results hold for the other two channels. 

The next set of three figures shows the same information for the screened Ohno potential 



2M given by 



V{^,J) = ^^^^^kMB (4) 



with the screening length (A) set to 0.75 of a C — C bond length(= a) (this A is somewhat 
larger than the Thomas- Fermi screening length for A^Cqq [^)- The major differences from 



the Hubbard model are that U* is generally higher, and that the region oi U* < has been 



pushed to higher U/t, as was predicted on the basis of second-order PT earlier p^j25|] . It 
is easy to see from Fig. ^ that the method is trustworthy only to about U/t ^3 — 4 now, 
especially for the Tig and Hg channels. Once again Fig. ^ shows the comparison between 
the different approximation methods. Adding a nearest neighbor interaction to the on-site 
Hubbard model leads to results very similar to Fig. ^. 

A simple physical picture of the reduction of U to U* is that the renormalized conduction 
level electrons interact weakly with the core primarily via the lowest Hg and Tig collective 
modes 0] (the second is actually a spin wave), which act very much like phonons in pro- 
ducing an effective interaction between the added electrons. While this picture produces 
qualitatively correct results, i.e., that the Ag pair state has the highest negative curva- 
ture, and that the Tig pair state always has positive curvature, the issue is complicated by 
nonzero overlaps between occupied and unoccupied levels in the interacting theory. We hope 
to present a full analysis in a future publication. 

To summarize, we have used a novel finite-size RG method to compute the effective 
Coulomb interactions at the FS for a single Cgo molecule. The method has an internal 
self-consistency check, and is reliable in the physically interesting range of bare interaction 
strengths, as opposed to other approximation schemes. The results show that while the 
effective interactions are always repulsive for physically relevant bare interaction strengths, 
there is a dramatic reduction in strength of the repulsion compared to its bare value. It is 
fascinating to speculate that for U/t > 6 in the Hubbard model two added electrons will 
attract each other on Cgo, although the method has reached the limit of its validity here. 
Clearly, this renormalization cannot be parametrized by the Anderson-Morel formula. 

A useful feature that we have not yet exploited is to reduce the size of the problem, 
especially for strong couplings. One could run the RG as long as the validity criterion 
allowed one to, and then exactly diagonalize the much smaller effective problem. 

The biggest deficiency of the above results is that since they include lattice screening 
phenomenologically at best (via Thomas- Fermi static screening |^), they cannot be applied 
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with much confidence to superconducting A^Cqo- We are currently working on a RG scheme 
which will include the effects of polarization and charge fluctuations of neighboring fuUerene 
molecules on U*, which should go some way towards addressing the lattice screening problem 
p3| . In the absence of such a calculation one can believe that lattice screening should only 
reduce U*, and that the viability of phonon mediated theories would be enhanced. In the 
Fig. 1^ we present a model calculation of the coupled electron-phonon system, where the 
effective interactions at the conduction level in the Hubbard model have been taken into 



account. We have used the phonon coupling constants of Ref [|12[- Despite the repulsive 
Coulomb interactions, the phonons succeed in binding the pairs. 

There are many directions in which this work could be taken. One could calculate the 



exciton energies |34] for neutral and charged Cgo, one could measure the spin structure of 



neutral Cgo, which has been the subject of some controversy [|35|, and one could hope to 
elucidate metallic screening on the lattice. 
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FIGURES 
FIG. 1. The Coulomb pseudopotential U* /t for the on-site Hubbard model as a function of 

bare coupling strength U /t for the Ag, Tig, and Hg channels. 

FIG. 2. The ratio of the two-loop to one- loop contributions at the final RG step as a function 
of U /t for the various channels for the on-site Hubbard model. Note that renormalized potential 
matrix elements are used in the loop calculations so that the result is reliable for much greater U /t 
than perturbation theory. 

FIG. 3. A comparison of different approximation methods for U* /t in the Ag channel. All the 
methods differ appreciably from our two-loop RG results long before the two-loop RG becomes 
unreliable. 

FIG. 4. The Coulomb pseudopotential U* /t for the Ohno potential with screening length 
A = 0.75a f« lA as a function of bare coupling strength U /t for the Ag, Tig, and Hg channels. 

FIG. 5. The ratio of the two-loop to one- loop contributions at the final RG step as a function 
of U /t for the various channels for the Ohno potential. 

FIG. 6. The same comparison as fig. 2 for the Ohno potential. Once again, all the methods 
differ appreciably from our two-loop RG results long before the two-loop RG becomes unreliable. 

FIG. 7. A model calculation of the binding energy of two electrons in eV , including the 
effective U* from the electronic calculation, and including up to two vibrons [|l8l. We have used 
t = 2 eV and the electron-phonon couplings from Ref. [0|. It is clear that even in the Ohno case 
the phonons succeed in binding the electrons. 
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